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Abstract 

In the presented review we discuss the problem of extraction the parameters of nuclear matter 
and charge distribution in stable and unstable isotopes. A substantial amount of information on 
the nuclear radii and other distribution parameters in light exotic nuclei has been obtained from 
experiments on intermediate-energy nucleus-nucleus collisions. The analyses of these experiments 
is usually performed with the help of the Glauber Theory. We review the main assumptions 
the Glauber Theory and show how this theoretical approach is used to calculate reaction and 
interaction cross sections. We also show how radii of nuclear matter can be obtained from the 
analysis of the experimental data on reaction cross sections. In the provided analysis reaction 
cross sections were calculated in optical and rigid approximations of the Glauber Theory as well 
as using explicit expressions. Numerical calculations of Glauber's explicit expressions for reaction 
cross sections were done using Monte Carlo technique. Recent results of the precise measurements 
of charge radii of light exotic nuclei which were done using the laser-spectroscopy technique as well 
as the method of the investigation of the nuclear matter distributions in proton elastic scattering 
experiments in inverse kinematics are also discussed. 
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I. INTRODUCTION 



The root-mean-square radii of nuclear matter (R^J as well as nuclear matter and charge 
distributions contain an important insight on nuclear potentials and nuclear wave functions. 

It seems that Elton p] was the first to carry out exact calculations of the influence of a 
finite size nucleus on the elastic scattering of high-energy electrons. He presented calculations 
of the ratios of differential cross sections for elastic eA scattering on extended and point-like 
nuclei, and showed that in some cases these ratios were very different from unity. 

The first detailed information about nuclear charge density distributions in stable nuclei 
comes from the experiments of R. Hofstadter. Many of his original papers together with 
other works concerning the discussed topic, including the paper of Elton, can be found in 



The main results of Hofstadter and results from other experiments are presented in his 
review j3]. Various models for charge density distributions p(r) were used to describe ex- 
perimental results. The list includes the point-like model p(r) = <5(r), uniform, Gaussian, 
exponential and some other distributions. First of all, it was observed that nuclei are char- 
acterized by a form factor S(q) 



which decrease rather fast with q 2 . Therefore, the point-like model can be excluded from 
further discussion. 

Secondly, nuclear matter is not distributed uniformly even in spherical nuclei. The radial 
dependence of the distribution is somewhat complicated function of the radial variable r, 
the distance from the center of the nucleus. Moreover, these distributions are different in 
different ranges of atomic weights A and charges Z. 

To describe the size of a nucleus the "effective" nuclear radius Ra was introduced. The 
experimental data show that for not very light stable nuclei its "effective" radius depends 
on atomic weight number A (see, for example, in [I]) 



The value of Rq slightly decreases when the atomic weight number A is increasing. 
The total inelastic hadron-nucleus cross section can be expressed in a following form: 



0- 




(1.1) 



R A = R A 1/3 , R = 1.2 + 1.4fm. 



(1.2) 



41 = *R\- 



(1.3) 
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The simplest expression for the total inelastic (or reaction) cross section, a AB , for interaction 
between nucleus A (with radius R A ) and nucleus B (with radius Rb) reads 

= ART + RT) 2 - (1.4) 

More accurate expression was suggested in [5] by Bradt-Peters 

C7« =7^(^/3 + 5V3 _ c)2j (L5) 

in which the probability of peripheral collision without nuclear interaction is taken into 
account. This expression is in reasonable agreement with the experimental data obtained in 
experiments on interactions of stable nuclei at high energy when R = 1.48 ± 0.03 fm and 
c = 1.32 ± 0.05, see |6j. However, it is clear that in the case of nuclei with halo the values 
of these parameters can be significantly different. 

Nuclear density distribution can be characterized by a set of its moments 

(rf) = J r 2n PA (r)d 3 r. (1.6) 

These moments are suitable when one compares different nuclear density distributions. All 
odd moments are zero due to obvious symmetry of density distributions. 

The zero momentum (n = 0) is equal to unity due to the normalization condition. The 
second momentum (n = 1) determines the root-mean-square nuclear radius, -R rms , is defined 
as 

R ms = , (r\) = J r' PA (r)d 3 r . (1.7) 

Generally speaking, the full set of all moments (r^ 1 ), n = 1, 2, ... determines the distribution 
p A (r) for every value of r. However, i? rms by itself is often used to compare nuclei described 
by different nuclear density distributions (see, for example, [7]). The higher moments are 
also discussed in the literature (see, for example, [8]). We will use notations R m and R c h 
when we consider separately the nuclear matter and nuclear charge r.m.s. radii. 

For stable nuclei the information on parameters of nuclear density distribution comes from 
the data on elastic scattering of fast particles on nuclear targets, [I]. The comparison of the 
data on electron and proton elastic scattering gives separate information about proton and 
neutron distributions in the nucleus, [TU] . The Glauber Theory, [TTHTB"] . is typically used 
to analyze experiments on interactions with nuclei at energies higher than several hundred 
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MeV. The obtained values of parameters for nuclear matter density and for charge density 
distributions are presented in [HI QUI EE] • 

In the case of unstable (radioactive) nuclei information comes from the experiments with 
radioactive beams, since the short lifetime often doesn't allow to make target using those 
elements. The unstable nuclei are produced through the projectile fragmentation of primary 
nuclear beam ( n B in [15], 18 in [16 J , 40 Ar in [17]) on the production target. The exper- 
imental production cross sections for different fragments produced in 40 Ar collisions with 
Be target can be found in [18]. The produced fragments are then separated by a magnetic 
analyzing system and rescattered on the nuclear reaction target (Be, C or Al). Again, the 
Glauber Theory is used to analyze these experiments at energies higher than several hun- 
dred MeV. However, the analytical calculation of all Glauber diagrams for nucleus-nucleus 
interactions is impossible and some approximate approaches can be used. 

In the presented review we mainly focus on the experimental results on nuclear den- 
sity distribution for unstable nuclei and on the problems of extracting parameters of their 
distributions. 

In Section [IT] we discuss the main assumptions of the Glauber Theory and its application 
to cross section calculations. 

The numerical results for the nuclear radii depend on the used approximation of the 
Glauber Theory. It is illustrated in Section |III| for the stable light nuclei, where the appli- 
cability of the Glauber theory and the nuclear matter distributions are well established. 

Historically, the first experimental study of unstable nuclei at high energies (790 
MeV/nucleon) was presented in [15l [19]. The interaction cross sections of He, Li and Be 
isotopes on Be, C, and Al targets were measured. It was found (independently on some 
analyses problems) that n Li has a radius much larger than other neighboring nuclei. This 
remarkable result suggests the existence of a long tail in the nuclear matter distribution in 
11 Li, i.e. halo. The obtained data stimulated the appearance of many experimental and 
theoretical papers in which various approaches were used to analyze experimental results 
obtained in the collisions of unstable nuclear beams with nuclear targets. The analyses of 



the experimental data and the main results are discussed in Sections IV 



In Section [V] we give a short review of the results for electrical charge distributions 
obtained laser-spectroscopy experiments. 

It is necessary to note that experimental measurements of the interaction cross section can 
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provide information on the only one parameter of the nuclear density distribution, or on the 
only one momentum of the distribution (for example, R rms )- Different type of experiments 
in which the differential cross sections for proton elastic scattering da/dt, where \t\ = q 2 , 
on unstable nuclei in inverse kinematics are measured, provide better information on the 
nuclear density distribution and several moments depending on the accuracy of experiment. 

In these experiments beam of unstable nuclei is scattered on the proton target and several 
experimental points at different values of squared transfer momentum q 2 can be obtained. 
That allows to determine separately, for example, the radii of nuclear core and halo. First 
results obtained using this method at projectile energies of about 700 MeV/nucleon were 
presented in [16]. Data analysis which provided information on core and halo radii was 



conducted in (20] (see Fig. 8). In Section VI we discuss results obtained in the experiments 
in inverse kinematics. 



Finally, a short summary and conclusion remarks are presented in Section VII 



Due to page limitation of this review we will not discuss the Coulomb contributions to 
the processes of light ion inelastic collisions. In the case of stable nuclei these contributions 
are small, about 3 ± 1% [2H 22] even for 12 C-Pb interaction cross section. They can be 
significantly more important for halo nuclei. For example the analyses of [21] estimates 
them to be 10 ± 3% for n Li-Al and 30 ± 12% for n Li-Pb collisions. 



II. GLAUBER THEORY FOR ELASTIC NUCLEUS-NUCLEUS SCATTERING 
A. Glauber Theory and its main assumptions 

Glauber Theory [TTHT3] describes interaction between high-energy particles and nuclei. 
The accuracy of Glauber Theory at intermediate energies was discussed in detail in [23]. The 
foundation of the Glauber Theory is the eikonal approximation for fast particle scattering 
in quantum mechanics. Description of the eikonal approximation can be found in [23], [22] 
or elsewhere. 

Let a fast particle (nucleon) with mass m, momentum k and kinetic energy T scatter on 
the nucleus A, which we treat for the moment as a collection of potential wells of the size a 
and of the depth Vq. Provided the following conditions 

ka > 1 , T/Vq > 1 (2.1) 
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are satisfied, the characteristic scattering angles are small. The phase shift XA(h) is given 
by the integral of the total nuclear potential (see detailed explanations in [21], [25J): 



oc 



Tfl 

XA(b) = --l dzV A (b,z), (2.2) 

K J — oo 



where b is the impact parameter (2-dimensional vector). Due to the conditions (2.1), the 
incident particle cannot interact more than once with a given target nucleon and the target 
nucleons have no time for interactions with each other during the scattering process. 

The next important assumption is that the nuclear potential V A is built up from spatially 
separated potentials which are corresponding to nucleons. In this case the phase shift on 
the nucleus, XA{b), can be represented as a sum of phase shifts for each nucleon- nucleon 
scattering, XN(bi)' 

A 

XA(b)=Y,XNN(bi) ■ (2.3) 
i=i 

In fact we do not have a good parametric characterization for the last assumption. Indeed, 
we saw that the distances between the neighboring nucleons are just of the order of the 
range of strong interactions. On the other hand, the small binding energy in comparison 
with nucleon or pion mass shows that, as a rule, nucleons within nuclei can be considered 
as quasi- free. 

B. Amplitude of elastic Nucleon-Nucleus scattering 



The phases Xa(P) and XNN(bi) in Eq. (2.3) are directly related to the elastic nucleon- 



nucleon scattering amplitudes on the nucleus and the isolated nucleon respectively: 

f$NiQ) = I d 2 bT(b)e^ b , (2.4) 



r NN (bi) = 1 - e*™^ = — j d 2 qe-^ f% N (q) , (2-5) 

r A (6;ri, ...,r A ) = 1 - e lXA(b;ri '-' rA) = 1 - e i ^t=^ (b ~ hi) , (2.6) 

where r\, ...,va are the positions of the nucleons and bi are their transverse coordinates. 
Because of condition Yji=i r i — 0, only A — 1 of the nucleon coordinates are independent. 
This leads to the so-called center-of-mass motion correction. 
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Let us now consider the elastic nucleon-nucleus scattering with the transition of the 
ground state of the target nucleus A to itself, 



ik 
2^ 



d 2 be i ^ b (A\r A (b;r 1 ,...r A )\A) . 



(2.7) 



Using Eq. (2.6), we obtain 



f na(q) = ^/ d 2 be* qb J d 3 r 1 ...d 3 r A -p(r 1 ,...,r A ) x 



{2.1 



i-n 



i=l 



2nik 



d 2 qie 



-iq 1 (b-b i ) 



/jviv(gi) 



Here 6, is the transverse coordinate of the i-th nucleon, and p(r\, ...,r A ) is the probability 
density distribution to find nucleons with coordinates (ri,r 2 , ■■■,r A ). 

Assuming that the internucleon interactions are rather week, we can neglect the possible 
internucleon correlations |1 54j . In this case probability density distribution p{r\, ...,r A ) can 
be reduced to the product of one-particle densities p(r^) 

A 

p(r 1 ,...,r A ) = J[p(r i ), J d 3 r iP { ri ) = 1 . (2.9) 



i=i 



Using assumption Eq. (2.9) and integrating over the positions of the nucleons we can get 
the following expression for F$ A (q) 



F e J A (l) 



ik 
2^ 



d 2 be iqb 



1 



1 



2nik 



d 2 qie-^ b f% N (qi)S(qi, 



(2.10) 



where S(qi) is the one-particle nucleus form factor determined by Eq. (1.1). 



C. Amplitude of elastic Nucleus-Nucleus scattering 

In the framework of the Glauber Theory the elastic scattering amplitude of nucleus A on 
nucleus B with momentum transfer q can be expressed as [261 EZ] 

FUl) = d%e^(A\(B\T AB (b ]ru ...r A] r[,...r> B )\B)\A) . (2.11) 

Here k is the incident momentum of one nucleon in A-nucleus in laboratory frame, and b 
is an impact parameter. This expression is written in the frame, where B-nucleus is a fixed 
target. 
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Let us assume, similarly to Eq. (2.3 ), that the phase shift in the nucleus-nucleus scattering, 
XabQ)), is equal to the sum of phase shifts for each nucleon-nucleon scattering, XNN{h)- 
Integrating over the longitudinal coordinates of all nucleons we obtain 

FU<1) = ^ J dV b [l - S AB (b)] , (2.12) 

where 

S AB {b) = (A\(B\i^[l-T NN (b + u i -s j )}\ \B)\A), (2.13) 
where Ui and s,- are the transverse coordinates of nucleons in the nuclei A and B, respectively. 



Contrary to the case of nucleon- nucleus interaction, integral in Eq. (2.13) cannot be 



evaluated analytically even with the assumption Eq. (2.9) that nuclear densities p{r\, ....ta) 
in both A and B nuclei are the normalized products of one-nucleon densities p(rj). 

To make the problem manageable, one can retain only part of all contributions in the 
expansion of the product in Eq. (2.13), that corresponding to the contributions character- 
ized by large combinatorial factors. The leading graphs correspond to the so-called optical 
approximation, [28, 29J, in which one sums up the contributions with no more than one scat- 
tering for each nucleon. In other words, only those products of amplitudes Tnn(P + — Sj) 
in Eq. (2.13) are taken into accounts which have different indices i, j. This approximation 
corresponds to the summation of diagrams shown in Figs, la, lb, lc,... In this approximation 
diagram which describes ra-fold interaction has a combinatorial factor A C n BC n . To avoid 
crowding of lines in Fig. 1 we have only shown the nucleon participants from the nucleus A 
(upper dots) and B (lower dots) with the links standing for interacting amplitudes, and we 
have not plotted the nucleon-spectators. 

In the optical approximation averaging and (B\...\B) of the product ILj[l ~~ 

^NN^b + Ui — Sj)] can be substituted with the averaging of the [1 — T^^ib + Ui — Sj)] 

S%(b) = II<4OT - IW& + TH - \B)\A) . (2.14) 

hi 

The amplitude of elastic nucleon-nucleon scattering can be described by the standard 
parameterization 

/^(g) = ^expf-i/3g 2 V (2.15) 



4vr r V 2' 

Here is the total NN cross section, is the slope parameter of the differential NN 
cross-section dependence on q 2 . We neglect the real part of f^^{q)- 




FIG. 1: Diagrams of the interaction of two nuclei, A and B, taken into account. 



When both nuclei are not very light, using the standard assumptions of the multiple 
scattering theory, scattering matrix Sab can be written in a following form 



where 



T opt {b) 



Sf B {b) 



tot 
°_N_N_ 

47T/3 



i A 



A 



T opt {b) 



d bid 6 2 Ta(6i)T b (6 2 ) exp 



exp[-T opt ( 

(b + h 



b,) 2 



2/3 



with Ta and Tb 



T A (b) = A dzpA (r = Vb 2 + ; 

J — oo ^ 

Neglecting the NN interaction range f3 in comparison to the nuclear radii, we have 

tot c 

T op t(b) = ^J d%T A (b - bjTzih) . 



(2.16) 

(2.17) 
(2.18) 

(2.19) 



In the diagram language Eq. (2.13) accounts for all possible intermediate states of nu- 



cleons between the interactions, as it is shown in Fig. 2a, while the optical approximation, 
Eq. (2.17), would correspond to the interactions with only one pole (nuclear ground state) 
in the both A and B intermediate states, [27] (see Fig. 2c). 

Unfortunately, numerical calculations in [30] (see also [31] for the case of collisions of very 
light nuclei) demonstrate that the optical approximation is not accurate enough even for the 
integrated cross sections. The difference with the data amounts for ~ 10 — 15 % in cr^ and 
it is even greater for differential cross sections, [30] . This disagreement can be explained by 
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FIG. 2: Two- fold interaction of two nuclei in the multiple scattering theory (a), in the rigid target 
approximation (b) and in the optical approximation (c). 



the fact that series with smaller combinatorial factors in Eq. (2.13) give significant global 



corrections to the optical approximation results. As a matter of fact, the terms of the series 
are alternating in sign, so, due to the cancelations of terms with opposite signs, the final 
sums of these series can have very different values. Thus, some classes of diagrams with 
non-leading combinatorial factors give significant contributions to the final total value. 

The rigid target (or rigid projectile) approximation, described in [32j|33], is more explicit 
than the optical approximation. It corresponds to averaging (B\...\B) inside the product in 
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Eq. (pll: 



S± B iP) = (A\ { H\(B\[1 - T NN (b + 1H- Sj)]\B) \ \A) = [T ? 



rg 



ry 



(2.20) 



where 



T rg (b) 



A 
x exp 



d 2 biTA(bi) exp 



a 



tot 

NN 



An/3 



d 2 b 2 T B (b 2 )x 



(6 + 61 



b,) 2 



2/3 



(2.21) 



Neglecting the NN interaction range in comparison to the nuclear radii T rg (b) can be sim- 
plified as 



T rg(b) = -^J d 2 b{T A (bi - b) exp 



tot 
a NN 



T B {h 



(2.22) 



This approximation corresponds to the sum of the diagrams in Figs, la, lb, lc, ... and the 
correction diagrams in Figs. Id, le, If, lg,... Diagrams in Figs. Id, le, If, lg,... represent the 
case when each nucleon from the nucleus A can interact several times, but all interacting nu- 
cleons from B are still different. Each correction diagram which describes n-fold interaction 
has a combinatorial factor smaller than the combinatorial factor leading diagrams (Fig. la, 
lb, lc, ...). Although, due to the obvious asymmetry in contributions of the two nuclei such 
approach can be theoretically justified in the limit A/B <C 1, i.e. for C-Pb collisions, this 
approximation can be used sometimes in the case of heavy ion collisions with equal atomic 
weights. 

It was suggested in [34J to present an expressions for the nucleus-nucleus scattering as a 
standard Glauber picture, where the amplitude of nucleon-target scattering is considered as 
an elementary particle and, after that, all Glauber rescatterings are taken in account. In [22] 
it was shown that such a picture is equivalent to the rigid target, or to the rigid projectile 
approach. 

Further corrections to the elastic amplitude (some of them are shown in Figs, lh, li, 
lj, lk,...) have been considered in [3"oTl3"§] . However, the results of such corrections are 
somewhat complicated for practical use. 

The possibility to obtain the Glauber Theory results without any simplification by the 



direct calculation of Eq. (2.13) using Monte Carlo simulation was first suggested in [HUE]. 
This method was used for numerical calculations in [351 B2HH] - The simplest algorithm 
considering the values of coordinates uniformly distributed in the interaction region cannot 
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be applied here because in most cases several coordinates have values corresponding to very- 
small nuclear density. 

The algorithm proposed by Metropolis et al in [45] allows to generate a set of coordinates 
which are distributed according to a pre-defined distribution. The Metropolis method (often 
reffered as Metropolis-Hastings [46] algorithm) is as follows: 

1. The initial coordinate Sj is randomly generated from the appropriate interval 

2. To obtain next coordinate, a shift Asj is randomly generated and then added to the 
initial coordinate 

3. New coordinate is accepted when the ratio r = p(sj + Asj)/p(sj) > 1 

4. If the ratio r < 1, new coordinate is accepted only when r > x, where x is a new 
random number x from [0, 1] interval. Otherwise new coordinate is not accepted 

Generated set of coordinates can be used to calculate average value of n«[l ~ ^NN^b + Ui — 
Sj)] and, Sab (b) . Some results of cross section calculations using described technique are 
discussed in Section [III 

D. Difference between reaction and interaction cross sections for light ion colli- 
sions 

The total inelastic (reaction) cross section for the collisions of nuclei A and B, <J AB , is 
equal to the difference of the total interaction cross section (J^ AB 

a% = ~ImF$ B (q = 0) = 2 J <fb[l - S AB (b)} (2.23) 

and integrated elastic scattering cross section cr A l B : 

<b = J ' d 2 b[l-S AB {b)f . (2.24) 

So, for the reaction cross section we obtain 

41 = <b - <b = J d 2 b[l - \S AB (b)\ 2 ] (2.25) 

The values of S A s(b) in these expressions can be calculated in one of the approaches described 
in Section ITTTl 
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As it was mentioned in Introduction, in the high energy light ion scattering experiments 
only interaction cross sections (J £1X6 measured. The difference between interaction and 
reaction cross sections is that the reaction cross sections include the cross sections of all 
processes except of the elastic scattering AB — > AB, whereas the interaction cross sections 
do not include the elastic scattering AB — > AB as well as the processes with a target nuclei 
excitation or disintegration AB — > AB* (B* ^ B), so a)[ B < 0"ab- The difference between 
a^B an d o"ab was estimated in [HI 09] to be not larger than a few percents of their values 
for a beam energy higher than several hundred MeV per nucleon. Usually this difference is 
neglected in the analyses of experimental data, see, for example, [3SJ 113 ED]- 

Lets define A' and B' as all excitation or disintegration states of A and B nucleons 
including their ground states, and A* and B* as all excitation or disintegration states ex- 
cluding their ground states. Therefore, difference between a A f B and a AB can be expressed 
as a^l - a^l = a(AB ->• AB*). 

The cross sections of the processes AB — > AB*, as well as the processes AB — > A*B*, 
where both nuclei can be exited or disintegrated, can be calculated within the Glauber 



Theory using the same assumptions as in Section II A 



Let us first calculate the cross section of the processes AB — > A*B*. The processes when 
an incident nucleus is excited [155] without changing its Z and N numbers are also very useful 
for our analysis of the difference between reaction and interaction cross sections, since they 
are not included into a^ B cross section. In the case of AB — > A'B', the amplitude of AB 
scattering together with the excitation or disintegration of one or both nuclei can be written 



in a form similar to Eq. (2.11) 

FAB-+A>B>(q) = ^J d 2 be i ^(A\(B\V AB (b;r 1 ,...r A ;r' 1 ,...r' B )\B')\A') . (2.26) 
Since all processes of AB — > A'B' transitions are the results of elastic NN scattering the 



operator Tab(&; •••?"a; r[, ■■■^b) has the same form as in Eq. (2.11 ). 

The difference between the total cross section a^ B and a(AB — > A'B') determines the 
cross section of secondary particle (pion) production, o v AB ! . Here we neglect the contributions 
of the processes where a pion can be produced in one NN interaction and can be absorbed 
by another nucleon or nucleon pair. 
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The cross section of all processes AB — > A' B' 

d&AB-^A'B' 1 



d 2 q 



k 2 



AB^A'B 



(?)P 



(2.27) 



A'B' 



can be calculated, following method described in [54]. 

Using the completeness condition for nuclear wave functions 

the cross section of AB — > A' B' can be written as 

*AB-+A>B> = /[I - 25^(6) + : 



(2.28) 



where 



/ AB (6) = (A\(B\ J]!! " r^(& + Ui - 8j )] \B)\A) . 



(2.29) 



(2.30) 



The cross section of the processes AB — > A*B*, i.e. processes with the excitation or 
disintegration of one or both nuclei without the elastic scattering channel, can be calculated 
as a(AB -> A*B*) = a(AB ->■ A'B') - o{AB -> AB). Therefore, cr(AB -> A*£*) can be 
expressed in the following form 



°ab^a*b* = J \Iab(P) - S 2 AB {b)]d 2 b . 



(2.31) 



Iab(V) can be simplified using parameterization of elastic scattering nucleon-nucleon am- 



plitude Eq. (2.15) and the optical approximation 



I%{b) = exp (6) 



(2.32) 



where 



T:; t {b) 



27T/3 



d%d% r A (6i)r fl (6 2 ) 



2/3 _ 



(2.33) 



When the radius of NN interaction can be neglected the expression in Eq. (2.33) can be 
simplified even more [27] 



ITb(P) 



exp 



a™ N d 2 bT A (b 1 )T B (b-b 1 ] 



(2.34) 
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where (T l f} N is the nucleon-nucleon inelastic cross section. 

As we mentioned above the difference between reaction and interaction cross sections 
is defined by the a(AB -> AB*). It is easy to see that a(AB AB*) = a(AB 
AB') — a(AB — > AB), if the projectile nucleus A has no excitation states. The amplitude 
of AB — > AB' processes has the form 

F A B^AB>(q) = ^f d%e^ b (A\(B\T AB (b ] r 1 ,...r A ;r' 1 ,...r' B )\B')\A) . (2.35) 

Again, using the completeness relation for the target B 

J2\AB')(AB'\ = l[\A)(A\5(s j -s f j ), (2.36) 



B' 



we obtain 



where 



*ab-+ab> = /[I - 2S AB (b) + J AB (b)}d 2 b , (2.37) 



J AB (b) = (A\(B\ I Y[[l - T NN (b + u t - Sj )] \ \A) x 



1,3 



X 



(A\ I]! 1 - T N N (b + u t - sj)] \B)\A) . (2.38) 



3 

h3' 



Using expression for the elastic cross section Eq. (2.24) we obtain following expression 
for a{AB -> AB*) 

o ab ^ab* = J [J AB {b) - S 2 AB (b)}d 2 b (2.39) 
In the optical approximation J AB can be expressed as 

JTb = exp(-2^(6)) , (2.40) 

where 

r;,(6) = / dXdXT^Ts^e-'^V 1 x 

xfl-^^/^We-^). (2.41) 

As usual we are using the parameterization of elastic scattering nucleon-nucleon amplitude 
Eq. (|2T5b. 
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If the radius of NN interaction can be neglected T* pt can be simplified 

TUb) = af N J d 2 bT A (h)T B (b - h) (l - ^T B (b - . (2.42) 

The difference between a AB and a AB for radioactive isotope 34 C1 and 12 C target was 
estimated experimentally in [UJ. It was found to be approximately 10 mb (i.e. about 1% of 
<?ab) an d in agreement with the calculations presented in [19]. The analyses in [55] results in 
larger difference between a AB and cr^B- ^ was f° un d that for experiments where projectiles 
with with A < 80 interacted with a carbon target the average difference between interaction 
and reaction cross sections was approximately 60 mb, i.e. about 4-6% of cr AB . However, 
provided in [55J results are not accurate (especially for the light nuclei, see in [56]), since the 
"black-sphere" model with rectangular distribution of nuclear matter was used for analyses 



E. One or several nucleon removal cross sections for stable nuclei 

Glauber Theory allows to calculate cross sections of one or several nucleons removal from 
the projectile nucleus A in the process of nucleus-nucleus collision. 

As we mentioned above, assuming that the internucleon interactions are rather week, we 
can neglect the possible internucleon correlations and express the probability density distri- 



bution p(ri, ta) as a product of one-particle density distributions p(r^) (see Eq. (2.9)). 
In this case the cross sections of one or several nucleons removal can be calculated using 
the AGK [571 158] cutting rules technique for all Glauber diagrams shown in Fig. 1 and to 
consider all possible intermediate states. As a result, the cross sections of removal of one, 
two, three, etc. nucleons (T AB , <J ab , o"ab> e ^ c - 111 the arbitrary AB collisions have the forms 



°% = A{a% - a^ B ) , (2.43) 

= A[A ^ A {-^B + 2^ " 2ai\ B ) , (2.44) 
(3) _ A(A - 1) (A - 2) (r) (r) (r) (r) 

a AB — 2j \ a AB 6a A-l,B + 6a A-2,B a A-3,B) • 

Here we took into account all Glauber diagrams shown in Fig. 1, including loop diagram 
Fig. Ik without any multidimensional integration. 
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An important point is that the cross sections cr^lj B , cr^*i 2 bi °"a-3 b i n Eqs. (2.43)-(2.45) 



cannot be taken from the experimental data. The cross sections in Eqs. (2.43)-(2.45) are 
results of the results of cancellation of the diagrams contributions (diagrams are shown 
in Fig. lb, le, lh, etc), which contain the radii of nuclei with weights A, A-l, A-2 in 



denominators. Therefore, cross sections <J^_ lB , u^^bi a A-3B should be calculated in the 

(r) 

framework of the Glauber Theory with the same nuclear radii, as cr\ B . 



The equations (2.43)-(2.45) are exact in the framework of the Glauber Theory. However, 
their application for the calculation of the physical processes, for example, to the cross 
sections of removal of several nucleons with the condition that the nuclear remnant stay to 
be bound needs in some additional physical assumptions. 

Really, we assume that in the process of NN scattering with transfer momentum of the 
order of that in free NN interaction both nucleons will be removed from the nuclei but the 
nuclear remnants will stay bound. It seems to be true for the collisions of light nuclei with 
an accuracy about 80% [60] 



The case of the fragmentation of unstable nuclei is discussed below in Section VC 



III. MONTE CARLO SIMULATION OF REACTION NUCLEUS-NUCLEUS 
CROSS SECTIONS 

A. Parameters of NN amplitude 

(r) 

In this section we will discuss numerical method to calculate reaction cross section a\ B . 
Below we present the results of the analysis the data obtained in the at high energy nucleus- 
nucleus collision experiments at energies about 800 MeV per nucleon. The parameters of 
NN elastic scattering amplitude at these energies were analyzed in [HI] . At energy 800 MeV 
per nucleon we obtained 

/3 2 =5 GeV~ 2 = 0.2 fermi 2 
a™ = 47 mb; oJ£ = 38 mb (3.1) 
a e J p = 25.4 mb; a e J n = 26 mb 

For the comparison and possible interpolation we present the same quantities at energy 
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1 GeV per nucleon: 

/3 2 = 5.5 GeV" 2 = 0.22 fermi 2 

a™ = 47.5 mb; a™ = 38.2 mb (3.2) 

a%, = 24.65 mb; a e p l n = 23.956 mb 

The ratio of the real to imaginary part of elastic scattering amplitude for p 4 He was 
obtained to be equal 0.06 ± 0.06 at energy 700 MeV, Since the term proportional to 
the second dergee of this ratio contributes to light ion cross sections, so it can be neglected 
in the further calculations. 



B. Extraction of R m values in various theoretical approaches 

As we mentioned in Introduction, for stable nuclei the information on parameters of 
nuclear density distribution comes from the data on elastic scattering of fast particles on 
nuclear targets. 

The total inelastic (reaction) cross section for the collisions of nuclei A and B, a^ B is 



shown in Eq. (2.23): 



= Jd 2 b[l-\S AB (b)f 



For the numerical calculations of Sab it is necessary to use an expression for the nuclear 
matter density distributions in the colliding nuclei. The most detailed information about 
these distributions comes from the data on electron or hadron differential elastic scattering 
cross sections on nuclear targets. 

The nucleon density in the light (A< 20) nuclei can be described by a harmonic oscillator 
(HO) density distribution [41147]: 

A/2 - 2 /r\ 2 \ / r 2 ' 

i— bJ exp -p 



w(r)=Pi(l + -4— ( T ) )exp(-- , (3.3) 



where A is the nucleus size parameter and p\ is the normalization constant. The nuclear 
density distributions in not very light nuclei can be reasonably described by Woods-Saxon 
expression 

Ur) = — VcV (3 " 4) 
l + e X p^j 



18 



Here p\ is the normalization constant, c is a parameter measuring the nuclear size, and a is 
related to the diffuseness of the nuclear surface. The parameter c shows the value of r at 
which p(r) decreases by a factor 2 compared to p(r = 0), p(r = c) = |p(r = 0). The value 
of a determines the distance r = 4a In 3 ~ 4.4a at which p(r) decreases from 0.9p(r = 0) to 
0.1p(r = 0). 

It was mentioned above that from the experimental data on a^ B it is possible to determine 



only one parameter of nuclear matter distribution, let us say R m , Eq. (1.7). It is necessary 
to note that the value of R m is somewhat smaller than the effective radius Ra — 1.2 A 1 ^ 3 fm. 
For example, in the case of uniform nuclear density with radius Ra, Ra — \J\( r A)- 

The reaction cross section of 12 C- 12 C interaction as a function of the R m of the 12 C nucleus 
was calculated in different approximations of the Glauber Theory for various nucleon density 
distributions. The results of the calculations are shown in Figure 3. 

In these calculations, which correspond to energies of 800-1000 MeV per projectile nu- 
cleon, the total NN cross section value, averaged over pp and pn interactions, = 43 
mb was used. The nuclear matter distribution parameters A, for harmonic oscillator (HO) 



density in Eq. (3.3), and c, for the Woods-Saxon density in Eq. (3.4), have been fitted to 



obtain the required R m , whereas the parameter a in Eq. (3.4) has been fixed to the value 
a = 0.54 fm. 

As the first step we have reproduced the result from [47] for a^ c _ 12c in the optical 
approximation without NN interaction range, using the harmonic oscillator (HO) density. 
Dependence of the reaction cross section, <r^2 C _\2 C , on the R m is shown by curve 2 of Fig. 3. 
It is in a good agreement with the result from |47], which is shown by the marker. 

Curve 1 in Fig. 3 represents the calculations of the dependence of the reaction cross section 
( y < \2 C _\i c on R m , which was done in the optical approximation with the Woods-Saxon density 
distribution. By comparing curves 1 and 2 of Fig. 3 one can see that equal R m with different 
assumptions about the nuclear density distribution result in different reaction cross sections. 
In other words this equivalently means that the same experimental reaction (or interaction) 
cross section with different assumptions about nuclear density distribution result in different 
R m , e.g. the assumption of the Woods-Saxon density distribution in 12 C nucleus leads to a 
smaller value of R m than the one obtained with harmonic oscillator density distribution. 

The corresponding results obtained in the rigid target approximation are shown by curves 
4 and 3 in Fig. 3. Here again the assumption of Woods-Saxon density distribution results in 
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2.10 2.20 2.30 2.40 2.50 2.60 2.70 
R ms > fm 



FIG. 3: The cross section of the reaction 12 C— C as a function of the R rms of the C nucleus 
calculated in different approximations with two different nucleon densities: 1) the optical approx- 
imation without range of NN interaction with Woods-Saxon density distribution; 2) the optical 
approximation without range of NN interaction with HO-potential density distribution; 3) the rigid 
target approximation without range of NN interaction with Woods-Saxon density distribution; 4) 
the rigid target approximation without range of NN interaction with HO-potential density distri- 
bution; 5) Glauber calculation with Woods-Saxon density distribution; 6) Glauber calculation with 
HO-potential density distribution. 

a smaller value of R m than the one calculated with the HO density distribution. The rigid 
target approximation contains additional diagrams Figs. Id, le, If, lg,..., which increase the 
shadow effects. That explains why both curves 3 and 4 lie below curves 1 and 2. 

The curves 5 and 6 in Fig. 3 show the results of the calculation of all the diagrams of the 
Glauber Theory by using Monte Carlo method and accounting for the finite range of NN 
interaction. Curves 5 and 6 were calculated with the Woods-Saxon density distribution and 
the HO density respectively. Here, new shadow corrections of the type shown in Figs, lh, li, 
lj, lk,..., appear by comparison to the rigid target approximation. As a result, the calculated 
reaction cross section becomes now smaller at the same value of R m . 

Clearly, the optical-limit approximation overestimates the calculated nucleus-nucleus re- 
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action cross sections. The difference between the reaction cross section calculated in the 
optical-limit approximation and that calculated with the help of the exact Glauber formula 
becomes even larger in the case of halo nuclei [32] • Therefore, to extract more accurate 
information on the nuclear sizes from the reaction cross sections, it is important to to per- 
form calculations with the exact Glauber formula. It should be admitted however that the 
calculations of the nucleus-nucleus reaction cross sections by the exact Glauber formula re- 
quire significantly more time compare to calculations in optical-limit approximation. The 
reaction cross sections calculations in the rigid-target approximation, being very simple, give 
noticeably more accurate results compare to the optical-limit approximation. 

The values of the interaction cross sections ail c _i2 C presented in [47] are 856 ± 9 mb 
and 853 ± 6 mb at energies 790 MeV and 950 MeV per nucleon, respectively. The older 
experimental measurement in [51] gives a value av2 C _i2 C = 939 ± 49 mb at energy 870 MeV 
per nucleon. The total 12 C- 12 C cross section was measured to be 1254 ± 54 mb at the 
same energy, whereas the Glauber Theory with R m taken from data in [47] predicts a value 
^c_i2 C = 1405 mb. 

In Table [T] we present the values of R m extracted from the measurements [I7J H7] of 
interaction cross section [T56] of stable projectile nuclei with l2 C target at energies 800-1000 
MeV per nucleon. 

The R m values were calculated assuming that the nuclear matter density distribution 



can be described by the Woods-Saxon expression Eq. (3.4). To get the dependence of the 
reaction cross section a^ B on the R m , we varied parameter c of the density distribution and 
kept parameter a as a constant at a = 0.54 fm. The R m values were extracted from the 
agreement of the calculated value with the experimental values of cr^ B . 

In Table [I] one can see that our calculations in the optical approximation with the Woods- 
Saxon density distribution and neglecting the NN interaction range result in the slightly 
smaller values of R m (0.05—0.1 fm) than those obtained in [UJ S7J. R m values are getting 
even smaller when calculated with a finite range of NN interaction. In the case of the Glauber 
Theory with the Woods-Saxon density distribution shadow corrections lead to larger values 
of R m than in the other calculations. 
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TABLE I: The values of R m in fm extracted from the measurements of interaction cross section in 
collisions of projectile nuclear beam with 12 C target at energies 800-1000 MeV per nucleon. Data 
on measured interaction cross section were taken from |17|, 14 7] . 

Nucleus Without NN range With NN range Glauber Theory 





oh, pang 


Optical 


Optical 


Rigid target 




C 12 


2.31 ±0.02 


2.25 ±0.01 


2.09 ±0.01 


2.18 ±0.01 


2.49 ±0.01 


N 14 


2.47 ±0.03 


2.42 ±0.03 


2.23 ±0.03 


2.35 ±0.04 


2.64 ±0.03 


16 


2.54 ±0.02 


2.48 ±0.02 


2.29 ±0.02 


2.41 ±0.03 


2.69 ±0.02 


F 19 


2.61 ±0.07 


2.55 ±0.08 


2.34 ±0.08 


2.44 ±0.09 


2.75 ±0.07 


Ne 20 


2.87 ±0.03 


2.84 ±0.04 


2.63 ±0.03 


2.75 ±0.04 


2.99 ±0.03 


Na 23 


2.83 ±0.03 


2.73 ±0.04 


2.52 ±0.04 


2.62 ±0.04 


2.91 ±0.03 


Mg 24 


2.79 ±0.15 


2.65 ±0.23 


2.44 ±0.22 


2.53 ±0.24 


2.85 ±0.20 


CI 35 


3.045 ±0.037 2.92 ±0.04 


2.68 ±0.04 


2.76 ± 0.04 


3.08 ±0.04 


Ar 40 


3.282 ±0.036 3.16 ±0.04 


2.90 ±0.03 


2.98 ± 0.04 


3.30 ±0.03 



C. Comparison of the calculated radii of matter distribution, R m , and radii of 
charge distribution, R c h 

Let us compare the values for R m obtained from nucleus-nucleus collisions and presented 
in Table [I] with the published results. It is known from [91 [10] that radii of proton and 
neutron distributions in nuclei with Z ~ A/ 2 are practically equal, so we can compare 
calculated radii for nuclear matter, R m , with electrical charge radii R c h presented in [1"4"] . 

It is necessary to make a distinction between the distributions of the centers of nucleons 
Pa(j) and the folded distributions Pa{t), where density pa(j) is convoluted with the matter 
or charge density in the nucleon, pjv(r): 

PA(r) = J p A (r - ri)piv(ri)<fVi . (3.5) 

We have to deal with pa(^) and with p^(r) when we calculate R m with and without ac- 
counting for the range of NN interaction, respectively. The rms radii of p^(r) and Pa{t), 
R m and R m , are different and the following relation between them was used in [9]: 

R^ = R^ ± (0.82 fm) 2 . (3.6) 
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In Fig. 4 we compare the result of R m calculations in the optical approximation (trian- 
gles) and in the rigid target approximation (squares) with values extracted from the 
electron-nucleus scattering experiment, [2]. Calculations of R m in both approximations 
were done with zero range of NN interaction (i.e. the folded distribution Pa{t) was used) 
and Woods-Saxon density distribution. Obtained results are systematically smaller than the 
data presented in [13]. This supports our point of view that in the optical approximation 
and in the rigid target approximation the effects of nuclear shadowing are too small, and, 
thus, one obtains agreement to the experimental nucleus-nucleus cross section with a smaller 
value of R m . 
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FIG. 4: The values of R c h extracted from electron- nucleus scattering (filled circles) and the values 
of R m obtained from nucleus-nucleus collisions in the optical approximation (triangles), and in the 
rigid target approximation (squares), both with zero range of NN interaction and with Woods- 
Saxon density distribution. 



The same calculations of R m were done in the framework of the Glauber Theory. The 
results of these calculations are presented in Fig. 5. In this case it is impossible to provide 
calculations with zero range of NN interaction because the contributions of diagrams with 
loops (see for example, Fig. Ik) present this range in the denominator. 
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The R m values calculated with distribution pAij) together with the R ch values extracted 
from electron-nucleus scattering experiment are presented in Fig. 5. They are in slightly 
better agreement to electron-nucleus scattering data than in the case shown in Fig. 4. 



• Electron scattering experimental data 
□ Glauber Theory with zero range of 

NN interaction • 

\ 

i 
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FIG. 5: The values of R c h extracted from electron-nucleus scattering (filled circles) and the values 
of R m obtained from nucleus-nucleus collisions in the Glauber Theory (open square) with Woods- 
Saxon density distribution and with finite range of NN interaction. 
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However, as we discussed before, the electron scattering data are related to the folded 
distributions pa{t) |157j . In order to make a more reasonable comparison, we calculate the 
needed R m in the case of the Glauber Theory as R m + A, where A was calculated as a 
difference between the R m obtained with the distribution of the nucleon centers and R m 



obtained with the nuclear matter density distribution, Eq. (3.4). Both root-mean-square 
radii were calculated in the optical approximation. 



IV. UNSTABLE NUCLEI, HALO AND SKIN 
A. Discovery of neutron halo 

The new era in nuclear physics was started in 1985, [T5l [19] . when the unstable nuclei 
with very large interaction cross sections were discovered. As the example of these sensation 
results we present in Table [n] the part of tables taken from [TS1 LH] on the interaction cross 
section of helium and lithium isotopes. 

It is evident that when the interaction cross sections for lithium isotopes with atomic 
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TABLE II: Interaction cross sections 07 of projectile nucleus on C target at 790 MeV/nucleon. 



null flf^l 1 Q 


/T r ( TY1 l \ 1 
u 1 I lllU J 




OOU zt 


4u p 
nc 


UUO HZ 




799 _l c 


8u p 


01 7 _i_ a 
ol 1 I U 


6 Li 


688 ± 10 


7 Li 


736 ±6 


8 Li 


768 ±9 


9 Li 


796 ±6 


n Li 


1040 ± 60 



mass in the range from 6 to 9 show some regular behavior, the cross section of n Li shows 
somewhat significant increase in its magnitude. Similar behavior was observed for helium 



isotopes. The simplest expression for the interaction cross section Eq. (1.4) together with 
the assumption that all these nuclei are spherical was used to analyze experimental data, 
|64j . The results taken from [12] for the radii of helium and lithium isotopes are presented 
in Fig. 6. The significant increase of n Li radius in comparison with nearest nuclei is evident 
and remains valid even after accounting for the criticism presented in [HI]. The classical 
behavior of nuclear radii A-dependences shown in Eq. (1.2) is evidently violated. 

The Hartree-Fock varionatial calculations with Skyrme potential for the structures of 
light nuclides are presented in [65J. It was shown that calculated interaction cross sections 
reproduce experimentally measured ones for all considered light nuclei except of 11 Li, where 
the calculated cross section is considerably smaller than the experimental one |158j . 

The experimental results obtained in [l~5l [19] and surprisingly large value for the n Li 
radius were interpreted in [67J as an evidence of a neutron halo existence appearing as a 
result of very low neutron binding energy of the valence neutrons. 

The nucleus u Li can be considered as a system of 9 Li core and two neutrons located in 
the halo - long tail of nuclear density distribution. In [HE] it was assumed that the density 
of n Li nucleus can be written as 

P^Li(r) = p c (r) + p h (r) , (4.1) 
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FIG. 6: Matter rms radii Rm . Lines connected isotopes are only guides for the eye. Differences 
in radii are seen for isobars with A = 6, 8 and 9. The n Li isotope has a much larger radius than 
other nuclei. 

where indices c and h are related to core and halo contributions. Using this assumption 
following values for R m values were obtained: R m ( 11 Li) = 3.1 ± 0.3 fm, R m (c) = 2.5 ±0.1 
fm and R m (h) = 4.8 ± 0.8 fm. 

6 He nucleus also can be considered as 4 He core and two neutrons located in the halo. 
Evidently, every neutron has a probability to belong to a core or a halo. 

Significant contribution of three-body forces between the core and two neutrons is a 
possible explanation of the anomaly large r.m.s. radius of n Li. An assumption that the n Li 
ground state should be considered as a three-body 9 Li+n+n system was confirmed later in 
many papers, see, for example, the analyses in [69], [TO] . An estimation of n Li wave function 
is presented in [70J. 

It was pointed out (see, for example, in [71]) that three particles interacting via short- 
range two-body interaction can form a variety of structures. The qualitative picture for the 
case of (A+n+n) system (core nuclide A and two neutrons) is shown in Fig. 7 (picture is 
taken from [80J). In such systems a so-called Borromean states can exist. Borromean state 
can be defined as a bound three-body system in which none of the two-body subsystems form 
a bound state. It is necessary to distinguish [72J Borromean states from Efimov's systems 
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FIG. 7: Schematic classification of three-body states of the system (A+n+n) as a function of the 
strengths of the two-body potentials Vau and V nn . The curve separates the regions where the 
three-body system is either bound or unbound, and the dashed lines separate bound and unbound 
two-body systems. The Borromean region is shown as dashed area. 

Both 6 He and n Li nuclei are considered as the excellent examples of a Borromean sys- 
tem confirmed experimentally. Both nuclei are stable in relation to the strong interaction, 
whereas nuclides 5 He and 10 Li as well as the dineutron nn state are the unbound systems. 

It is interesting to note that stable nuclei can have the Borromean excited states. As an 
example [71], let us consider the Hoyle resonance in 12 C (0 + state with energy 7.65 MeV). 12 C 
nucleus can be considered a cluster of three a-particles and every two-body system here is 
unbound 8 Be state. This Hoyle resonance plays very important role in stellar nucleosynthesis 
|75ti77] and, in some sense, the existence of our Universe is based on its existence. 

In addition to the Borromean states, the three-body system can exist in the configu- 
rations, first one where one two-body system is bound and two are unbound ("Tango" 
configuration, [78]) and another where two two-body system are bound and one is unbound 
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("Samba" configuration, [79]). It is possible, that 20 C nuclide is an example of "Samba" 
system, which is composed of two-neutron halo and a 18 C core. The neutron and 18 C forms 
the weakly bound state 19 C. 

B. Neutron (proton) halo and skin 

After the discovery of several unstable nuclei with unusual nuclear matter distribution, 
the existence of neutron or proton halo [151 IH] and neutron or proton skin [83J became 
clear. 

Usually, halo is considered as a long low density tail in the nuclear matter distribution 
[82] , whereas skin means a significant difference between R m values for protons and neutrons. 
The difference between these two scenarios depends on the values of some parameters which 
can be determined in theoretical Hartree-Fock calculations. Sometimes same nuclei are 
considered by different authors to have halo or skin. For example, in [83] authors described 
nuclides 6 He and 8 He as ones with nuclear skin, meanwhile in [84] same nuclei are described 
as nucleon with two-nucleons and four-nucleons halos. Very clear outlook of this problem 
was presented in [85j [86] . 

Firstly, let us consider a nucleus with given values of Z and N which nuclear density 
distribution exhibits the presence of neutron skin. Usually the ratio of neutron to proton 
density distributions, p n (r) / p p (r), is about N/Z in the interior region of a nucleus. This ratio 
can somewhat exceed N/Z near the nuclear surface, which can be explained by fluctuations 
due to shell effects, etc. The following criteria were introduced in [85] to define the neutron 
skin: 

1. in the neutron skin 



It was recommended in [85] to use £i = 4. 



2. contrary to the case of neutron halo, a neutron skin should contain a significant number 
of neutrons. 



Pn(r)/p P (r) > fi for r ~ R A . 



(4.2) 




Pn 



r ~ R A )/Pn(r = 0) > & . 



(4.3) 



Again, it was suggested in [HS] that the value of £ 2 should be 1/100 



28 



3. the difference between neutron and proton radii should be large enough, i.e. 



5R = R n — R„ > £3 



(4.4) 



where £3 = 1 fm [85] . 



The existence of nuclei with proton skin is also under discussions. The proton skin can 
be defined in the same manner as it was done for the neutron skin in Eqs.(4.2)-(4.4). 8 B 



and Al are amongst the possible nuclei with proton halo. In the case of Al the value 5R 



defined as in Eq. (4.4) is equal to -0.84 fm [83]. Generally, proton halos are equally possible 
as a neutron halos but less pronounced because of the strong influence of the Coulomb 
barrier. 

Definitions of halo and skin presented in |85j suggest that the main difference between 
nuclear skin and nuclear halo is the second condition, see Eq. (4.3). The skin is a phenomenon 
which involves rather large number of nucleons, whereas halo is caused by one or two nucleons 
extremely loosely bound with the core, that results in an abnormal slope of nuclear density 
distribution tail. 

The general outlook on light nuclei with halo or skin is presented in Fig. 8 taken from 

ESI. 
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FIG. 8: The light end of the chart of nuclides. Some of the drip line nuclides found to exhibit new 
phenomena, such as halos and skin. 
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It seems that the skin/halo effects decrease with increase of atomic weight. For example, 
the cross section or interaction radii measured in the collisions of Na isotopes on a carbon 
target [90] in a wide range of atomic weights (A=20-23, 25-32) increase with the number of 
neutrons rather regularly, contrary to the behavior of He and Li isotopes shown in Fig. 6. 
The neutron skin in neutron-rich Na isotopes is only slightly larger in comparison with skin 
in Ca isotopes, [ID] . 

The theoretical analysis provided in [80] results in the necessary, but not sufficient condi- 
tion for occurrence of two- or three-body halo systems (i.e., say, one- neutron or two-neutron 
halos) : 



where Sn is one-nucleon or two-nucleon separation energy and A is atomic weight. 

Detailed discussion of neutron halo nuclei can be found in [87]. The theory of nuclear 
halos can be found in [9_TJ [92] . 

The so-called few-body approach of the Glauber model was suggested in [9"3T[9l)] for the 
analyses of reactions for halo nuclei. The nucleus with halo, say 6 He or 11 Li, was considered 
as the systems of core (c), 4 He or 9 Li, and two valence (v) neutrons. The nuclear density 
distribution of the halo projectile (P) nuclide can be written as 



where the relative motion of core and valence neutrons is accounted for. If the interactions 
of core and valence nucleons are calculated in the optical limit of the Glauber Theory, 
the few-body approach theoretical expressions are different, |96j, from the standard optical 
limit, so the numerical results are also different. However, in the complete Glauber Theory 
the results of calculations in the few-body approach for 6 He- 12 C and 11 Li- 12 C interactions 
practically coincide up to numerically small contributions coming from nucleon correlations 
in the few-body approach [35J with the results of the standard Glauber expression. 

C. Nucleon separation energies and fragmentation cross sections for unstable nu- 
clei 

Let us consider an unstable nucleus as a core and one nucleon distributed at somewhat 
long distance from the core (in units of the range of the strong interaction). The probability 



S N A 2/3 < (2 + 4) MeV , 



(4.5) 



Pp{r) = p c (r)+p v {r) , 



(4.6) 
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p(r) to find such nucleoli at distant r from the centre of a system is defined in quantum 
mechanics [21], [2H] by the square modulus of the correspondent wave function 

p{r) ~ -e- 2 * , (4.7) 
r 

where 

fc = \/2me . (4.8) 

Here e is the bound energy and m is a nucleon mass. 

Evidently, the nuclides with halo should have a small separation energy of a neutron or 
a proton, S n or 5* p , in comparison with stable nuclei, where this energy is about 6-8 MeV 

The nucleus of n Li can be used as a clear illustration of this phenomenon. Usually n Li is 
considered as a core 9 Li with two neutron halo (long distance tail in the n Li wave function). 
The separation energy of two neutrons from 11 Li nuclide defined as 

S 2n ( n Li) = AE[( 9 Li + 2n) - ( n Li)] (4.9) 

is extraordinary small 

5 2n ( n Li) = 0.34 ± 0.04 MeV [HE] 
S 2 „( n Li) = 0.295 ± 0.035 MeV [US] 
S 2n ( n Li) = 369.15 ± 0.65 keV [TM] . 

The main candidates for being nuclei with neutron or proton halo together with their one- 



neutron and one-proton separation energies, and mean life times are presented in Table III 
In some cases, for example, for 6 He and n Li nuclides, the one-neutron and two-neutron 
separation energies coincide because the nuclides 5 He and 10 Li are unstable and the second 
neutron is emitted during a nuclear time-scale. 

Generally, the halo nuclei have one-nucleon or two-nucleons separation energies less or 
approximately 1 MeV. However, there are some exceptions. Amongst these exceptions are 
helium isotopes, which neutron separation energies are 1.9 MeV for 6 He and 2.6 MeV for 
8 He. That permits us to say that these nuclides have neutron skins. 

At the same time the nuclide 14 B has one-neutron separation energy S n = 0.97 MeV, but 
it is not considered as a halo state. Possibly, the reason is that in beta-decay of 14 B the modes 



(3 n and (3 nn are not observed, contrary to the cases of nuclides presented in Table |HI 
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TABLE III: Possible halo states. 



Nuclide 


S n (keV) 


S p (keV) 


r 


Configuration 


b He 


1867 ± 50 


26520 ± 95 


807 ± 15 ms 


4 He + 2n 


8 He 


2574 ± 18 




119 ± 12 ms 


4 He + 4n 


i i T . 

n Li 


325 ± 31 


15303 ±31 


8.59 ±0.14 ms 


9 Li + 2n 


n Be 


504 ± 6 


20165 ± 16 


13.8 ±0.8 s 


10 Be + n 


14 Be 


1850 ±51 




4.35 ±0.17 ms 


12 Be + 2n 


8 B 


13020 ± 70 


137 ±1 


770 ± 3 ms 


7 Be + p 


17 B 


1430 ± 150 




5.08 ± 0.95 ms 


15 B + 2n 


19 C 






49 ± 4 ms 


18 C + n 


17 F 


16800 ± 8 


600.3 ±0.2 


64.49 ±0.16 s 


16 + p 



where such decays exist. Even the halo state of 8 B where the proton separation energies 
is smaller than 0.14 MeV is not considered as evident, [82 J. All these facts indicate that 
today there is no the formal scheme which would allow to determine halo and skin nuclides. 
On the other hand, the two-neutron halo structure of 6 He, 11 Li, and 14 Be is successfully 
reproduced by theoretical calculations [79] . 

Let us consider now the cross sections of fragmentation of the unstable nucleus with halo 
(A) into fragment F, which is the core of A, after collision with stable nucleus (B) from a 
target. The fragmentation process can be expressed as AB — > F+X, where X consists of all 
possible states of halo nucleons together with possible production of secondaries. All states 
of target nucleus are also included in X. 

The nuclear density of a halo nucleus can be considered similarly to Eq. (4.6). Using 



expressions Eq. (2.16)-(2.19), the difference between reaction cross sections a^ B and af B 
can be written as 



a 



(r) 
AB 



a 



(r) 
FB 



d 2 b 



% J d 2 6T c (&-& 1 )T B (6i) 



X 



X 



/ d 2 6T fe (6-6i)T B (6i) 



(4.10) 



The obtained expression can be viewed on as probability of the process in which the core 
of nucleus A doesn't interact and nucleons from halo does not interact. The probability of 



the former is expressed as exponential factor in (4.10) and the probability of the latter is in 
square brackets. 
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So, it should correspond to the cross section of AB — > F+X fragmentation cross section 

a(AB ^ F + X) = , (4.11) 

in accordance with the result presented in [19] for the case when the beam nucleus is n Li 
and the fragment nucleus is 9 Li. 



In Table IV we compare the fragmentation cross sections with the differences of beam 



and fragment reaction cross sections. Since expression Eq.(4.11) for a(AB — > F + X) was 
obtained in the optical approximation without accounting for the shadow effects between 
core and halo distribution, it is hard to expect perfect match with experimental data. Nev- 
ertheless, one can see a reasonable agreement of the experimental data with calculations 
done using expression Eq. (4.11) for the cases of 14 Be — > 12 Be, n Be — > 10 Be, n Li — > 9 Li and 



6 He — )- 4 He fragmentation. 

The absence of agreement between experimental data and calculation done using 



Eq. (4.11) for 8 He — > 4 He and 8 He — )- 6 He fragmentation requires additional explanation. Af- 
ter interaction of one or several valence neutrons of 8 He with the target, two bound states, 
6 He and 4 He, can appear after the final state interactions. So, one can assume that the sum 
of two and four neutron separation cross sections of 8 He nuclide should be very close to the 
difference in interaction cross sections of 8 He and 4 He: 

(j( 8 He ^ 6 He) + (x( 8 He ^ 4 He) = 297 ± 19mb , 
a (/) ( 8 He) - a (/) ( 4 He) = 314 ± 8mb . 



V. MATTER AND CHARGE RMS RADII IN UNSTABLE NUCLEI 
A. Data at low energies 

There is a large volume of the experimental data on nucleus-nucleus interactions obtained 
at relatively low energies smaller than 300-400 MeV per nucleon. These data can be useful 
for measurements of the mass of nuclide, energy separation of a valence nucleon, etc. 

The data shows the qualitative difference in differential cross sections of elastic p 9 Li and 
p u Li scattering [97] . That provides additional arguments for an existence of neutron halo 
in n Li nucleus. However, independently on the quality of available data, the quantitative 
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TABLE IV: Production cross sections of the projectile fragments a(AB — > F ' + X) on 12 C target at 
790 MeV/nucleon and differences of beam and fragment reaction cross sections also on 12 C target. 



Beam Fragment 


a{AB -)• F + X) 


(r) M 
AB FB 






(mb) 


(mb) 


14 Be 


12 Be 


210 ± 10 a 


182 ± 71 a 


n Be 


10 Be 


169 ± 4 a 


129 ± 13 a 


n Li 


9 Li 


213 ±21 6 


260 ± 20 a 


8 He 


6 He 


202 ± 17 6 


95 ±9 a 


8 He 


4 He 


95 ±9 b 


314 ±8 a 


6 He 


4 He 


189 ± 14 6 


219 ±8 a 



"Data were taken from [102 . 
h Data were taken from [88] . 



interpretation of the data and the extraction of parameters of the nuclear density distribution 
(or R m ) for unstable nuclides are problematic. 

As it was mentioned in Section II A[ the corrections to the interaction and reaction cross 



sections calculated in the framework of Glauber Theory become significant in the low-energy 
range. For example, it was shown in [98] that taking into account a Fermi-motion of nucleons 
in projectile and target nuclei changes the effective energy of nucleon-nucleon interaction 
(NN). Since the cross section of AW interaction has a strong dependence on the effective 
energy, when at energies below 100 MeV per nucleon, the correction to the reaction cross 
section reaches 10% of its value. 

At smaller energies the modification of the fast proton trajectories by the Coulomb field 
leads to even more significant corrections. It was shown in [99] that at energy 30 MeV per 
nucleon this correction can change value of cr^ for 12 C - 40 Ca collision by approximately 
10%. Moreover, the mentioned above corrections are significant even for calculations of 
low-energy proton- nucleus inelastic cross sections, |100j . 

In order to extract parameters of nuclear density distribution the following semiemprirical 
parameterization of reaction (or interaction) cross sections was suggested in [101J: 



= nR 2 



A 1 / 3 A 1/3 

K /3 + A" + a Z w.s - C(E) 



4 /3 + a] /3 



B c 



E r 



(5.1) 
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Here Rq is radius of nucleon-nucleon interaction, A p and A t are the atomic weights of pro- 
jectile and target nuclei, a is an asymmetry parameter, C(E) is the energy-dependent trans- 
parency parameter which is constant at E > 100 MeV, and 

B - ZpZte2 f5 2) 

c "i.3(4 /3 +4 /3 ) 

is the Coulomb barrier. 



It was shown in 1 Of] that expression for the reaction (interaction) cross sections Eq. (5.1) 
is in good agreement with the experimental data obtained with stable nuclei. However, 
contrary to the general belief, this expression cannot be applied for the quantitative analyses 
of experimental data with nuclei near the drip lines, where nuclei with the halo or skin can 
exist, [T02] . 

First of all, the parameter C(E) which accounts for the diffuseness of nuclear surface 



(similarly to the parameter c in Eq. (1.5)) obtained from the analysis of interaction between 
stable nuclei cannot be used to describe data obtained in experiments with unstable nuclei. 
Clearly, the value for the diffuseness parameter for unstable nuclei with halo should be 
significantly higher compare to its value for stable nuclei. Secondly, semiempirical expression 



Eq. (5.1) suggests that the the value of Ro is the same for projectile and target nuclei, which 
is inconsistent with the experimental data obtained at low energy [105J. For example, the 



based on the Eq. (5.1) analysis would yield different values for the radius of the target 
nucleus 12 C when applied for the data obtained in 12 C- 12 C and n Li- 12 C collisions. As a 
result, anomalously large value for the radius of n Li nuclide becomes "distributed" between 
the radii of n Li and 12 C. 



For these reasons we believe that expression Eq. (5.1) cannot be used to analyze ex- 
perimental data with radioactive nuclei, and, hence, we will not present the values of R m 
obtained with the help of this expression in Table [V] in next section. 

Of course, all these criticism cannot be applied to the case of qualitative comparison of 
two projectile nuclides interacting with the same target. 
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B. Summary of results for nuclear radii extracted from interaction cross section 
data at high energies 



In this section we present the values for root-mean-square radii, R m , of nuclear matter 
distributions extracted from the experimental data on interaction cross sections obtained at 
somewhat high energies (higher than 600 MeV per nucleon). All data on rms radii presented 
in this section were obtained in optical approximation of the Glauber theory. Comparison 
with the data for nuclear radii presented in Table 1 shows that the difference between values 
for nuclear radius obtained using exact formulas of the Glauber Theory and values obtained 
in optical approximation is approximately the same for all nuclei. We do not present the 



values of interaction radii obtained by using Eq. (1.4) since these radii were obtained without 
taking into account the fact that the diffuseness region of the nuclear density distribution 



(parameter c in Eq. (1.5)) varies significantly when a halo present. 

The data |159] presented in [TH] were obtained using a shell-model harmonic oscillator 
function [1] as a nuclear matter density distributions. The results obtained after an im- 
provement of the data presented in [19] together with new results for Be and B isotopes are 
presented in |109] . In the later reference free values of NN cross sections were used and the 
R m values for point nucleon distributions are presented. Due to these reasons the R m values 
presented in [109] are significantly smaller in comparison with the results in [i~9] . 

In the case of 8 B nucleus there are some inconsistencies between somewhat small R m value 
obtained using the harmonic oscillator density distribution in [106J and very small separation 
energy of one proton |107j . S± p = 0.14 MeV, as well as with results of calculations, see in 
[106] . In [89] 8 B nucleus was considered as a possible candidate for having the proton halo, 
see Fig. 7. It is necessary to note that the R m value obtained in |106] for 8 Li nuclide was 
also smaller than the value presented in [19] , see Table |VJ 

The values of R m for carbon isotopes |110j as well as for 17 N, 17 F, and 17 Ne were obtained 
for point-like nucleons using harmonic oscillator nuclear matter density distribution |108j . 

The nuclear matter radii of A=20 isobars were analyzed in |lllj using Woods-Saxon 
nuclear density distributions. Results of this analysis are also presented in Table |V| 
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TABLE V: R m values of nuclear matter distributions extracted from the data on high energy 
(Ei a b > 600 MeV) nucleus- nucleus collisions. 



Nucleus R r mm (f™) 


Ref. 


Nucleus R rmm (fm) 


Ref. 


4tt q 

rie 


i 70 j- n 
i. i z ± u.uo 


[191 


14n„ 
Lie 


Q 1 £i 4_ n QQ 
O. iO ± U.OO 


[109] 


rie 


l.D f ± U.U4 


11091 


L> 


o a q -L n nQ 
Z.4o ± U.Uo 


|106j 


rie 


O Q7 _1_ n CiA 

z.o l ± U.U4 


nui 


8T3 

o 


o qq _i_ n n/i 
z.oo ± U.U4 




ne 


o /i q 4- n nQ 

Z.40 ± U.UO 


11091 


12 R 


o qq 4- n no 
z.oy ± u.uz 


[109] 


rie 


O Q 1 i n no 

z.oi ± U.Uo 


[T9l 


L> 


o a a 4_ n i o 
z.40 ± U.lz 


|109j 


8u„ 
rie 


o 4- n nQ 

Z.OZ ± U.Uo 


11091 


L> 


o /i/i 4- n nR 

z.44 ± U.UD 


[109 


Ll 


O £ /I 4_ n nQ 

Z.04 ± U.Uo 


[191 


L> 


O /I 7 4_ n 07 

Z.4 i ± U.Z / 


11091 


Ll 


O QO _1_ H HQ 

Z.oz ± U.Uo 


fT09l 


Ly 


o a q 4_ n no 
Z.4o ± U.Uz 


119 




Z.04 ± U.UO 


[191 


12p 


o qi 4- n no 
z.Ol ± U.Uz 


|110j 


Ll 


o qq _i_ n no 
z.oo ± U.Uz 


109 




O OQ 4_ n CiA 

z.zo ± U.U4 


[TTUl 


J_il 


o -i- n nQ 

Z.O ( ± U.UO 


[191 




o Qn 4- n r\T 

Z.ol) ± U.U / 


[110] 


1j1 


o Q7 ~L n no 
z.o ( ± U.Uz 


|109| 


16p 


o vn 4- n nQ 
z. / U ± U.Uo 


|110j 


8T i 


o Q7 -A- n no 

z.O I ± U.UZ 


[1061 


17p 


o 70 4- n nQ 

Z.lZx U.UO 


|110j 


9T i 


o kc\ -L n no 
z.OU ± U.Uz 


[191 


18p 


o qo 4- n n/i 
z.oz ± U.U4 


|110j 


9T i 


o qo j- n no 

Z.OZ ± U.UZ 


11091 


19p 


q i q 4- n n7 
o.io ± U.U / 


|110j 


Ht ; 
Lil 


Q Q£ _1_ n O/l 

o.oo ± U.z4 


[191 




o nQ 4_ n n £ 
z.yo ± U.UO 


[110] 


Ht ; 
Lil 


q i o _i_ n i fi 
o.lz ± U.iO 


109 


IN 


o a q 4_ n ncc 
z.4o ± U.UO 


108 


Lie 


n A 1 in no 
Z.41 ± U.Uo 


[191 


20t\t 
IN 


o 77 4_ n riA 
Z. 1 1 ± U.U4 


Ill 1| 


7 Be 


o qi _i_ n no 

_,.v.l_l. —1— L/ . VJ Zj 


11091 


20 q 

V J 


o «4 _i_ n nQ 


111 1| 


9 Be 


2.53 ±0.01 


[m 


17 F 


2.54 ±0.08 


[1081 


9 Be 


2.38 ±0.01 


[EH] 


20p 


2.75 ±0.03 


llil] 


10 Be 


2.43 ±0.02 


EE] 


17 Ne 


2.75 ±0.07 


(108] 


10 Be 


2.30 ±0.02 


[109] 


20 Ne 


2.84 ±0.03 


[111] 


n Be 


2.73 ±0.05 




20 Na 


2.69 ±0.03 


[111] 


12 Be 


2.59 ±0.06 


|109] 


20 Mg 


2.86 ±0.03 


|111| 
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C. Charge radii of light nuclei 



The classical way to measure the charge (or proton) radii of the nuclei is to analyze 
experimental data obtained in the eA elastic scattering experiments [3J. However, in the 
case of unstable nuclei this way is currently still under discussion |112l 1113] . 

The best information about charge radii of unstable nuclei comes from the laser spec- 
troscopy experiments. The energy levels (or the isotope shifts) of atomic electrons can be 
measured with very high accuracy, ~ 10~ 6 or better, |114[ II 15 j in relative units. Experi- 
mental results are in very good agreement with theoretical calculations. 

The isotope shifts between the atomic transitions in basic stable nuclei, 4 He, 7 Li and 9 Be 
and their unstable isotopes are analyzed. Some contributions to the calculated values of the 
shifts are proportional to the ratio of the finite charge radii of the stable nucleus and its 
unstable isotope |114[ I115j . That allows to extract the difference in the values of R c h in the 
basic stable nucleus and an isotope nuclide by using the results of theoretical calculations 



based on QED. The obtained results are presented in Table VI 



The radius of the point-like proton distribution, R p was estimated in |117] to be 

R p ( 6 He) = 1.912 ±0.018 fm . (5.3) 

One can see noticeable decrease of R c h in the intervals from 6 Li to 9 Li and from 7 Be to 
n Be and 10 Be and significant increase of charge radii for n Li and u Be. That confirms the 
unusual structure of the last nuclides. The 6 He, 8 He, n Li and n Be nuclei are supposed to 
consist of the 4 He, 4 He, 9 Li and 10 Be cores and halos of 2 neutrons, 4 neutrons, 2 neutrons 
and 1 neutron, respectively. However, the charge radii of 6 He, 8 He, n Li and n Be are larger 



than the corresponding radii of 4 He, 4 He, 9 Li and 10 Be (see data provided in Table VI). The 
larger charge radii of 6 He, 8 He, n Li and n Be can be explained by the core motion around 
the center-of-mass of these nuclei and partially by the possible nuclear core polarization 
[69| 1701 1125H128] . The contribution to the effective core size due to its motion around the 
nuclear center-of-mass can be determined from the data on the Coulomb nuclear dissociation 
[T2TM34] 

Some results on nuclear charged radii were extracted from the data on only interaction 
cross section. Since some additional assumptions on the nuclear structure are required to 
analyze experimental data, we will not present them in the article. 
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TABLE VI: R c h values of nuclear charge distributions obtained from laser spectroscopy data. 



Nucleus 


Rch (fm) 


Ref. 


Mean life 


4 He 


1.673 ± 0.001 




Stable 


6 He 


2.054 ±0.014 


1117] 


807 ms 


6 He 


2.068 ±0.011 


H18J 


- 


8 He 


1.93 ±0.03 


Iil8j 


119 ms 


6 Li 


2.51 ±0.06 


Ulfll 


Stable 


6 Li 


2.540 ±0.028 




Stable 


6 Li 


2.49 ± 0.04 


[120J 


- 


6 Li 


2.55 ±0.04 


imj 


- 


7 Li 


2.39 ±0.03 


imj 


Stable 


8 Li 


2.29 ±0.08 


[119J 


840 ms 


8 Li 


2.281 ±0.032 


H15J 


- 


9 Li 


2.22 ±0.09 


1119J 


178 ms 


9 Li 


2.185 ±0.033 


imj 


- 


9 Li 


2.217 ±0.035 


(1221 1123] 


- 


n Li 


2.426 ± 0.034 


1115J 


8.6 ms 


n Li 


2.467 ±0.037 


U23J 


- 


7 Be 


2.645 ±0.014 


[115J 


53.2 d 


7 Be 


2.647 ±0.017 


[124] 




9 Be 


2.519 ±0.012 


11241 


Stable 


10 Be 


2.357 ±0.018 




1.5 x 10 6 3 


n Be 


2.463 ±0.016 


[124] 


13.8 s 



VI. NUCLEAR MATTER RADII FROM PROTON-NUCLEUS ELASTIC SCAT- 
TERING AT INTERMEDIATE ENERGY 

Proton-nucleus elastic scattering at intermediate energy is an efficient means for study- 
ing nuclear matter distributions [10]. There is an obvious advantage of proton scattering 
experiments at intermediate energy as compared to similar experiments at low energy. As 
we discussed before, at intermediate energy the mechanism of proton-nucleus scattering is 
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somewhat simple and can be described in the framework of various modifications of multiple 
scattering theories, in particular the Glauber theory [13]. That allows one to connect rather 
accurately the measured differential scattering cross sections with the nuclear matter distri- 
butions under study. A number of experiments on proton-nucleus scattering were performed 
previously to study the matter distributions in stable nuclei. Such experiments can be also 
carried out to study matter distributions in short-lived unstable nuclei. In this case the 
experiments should be performed in inverse kinematics. 

The scattering of protons from the nuclear halo is confined to small scattering angles. 
Therefore, in order to study the spatial structure of halo nuclei, it is important to measure the 
differential cross sections for proton scattering at small momentum transfers. The analysis 
of the differential cross sections for proton scattering at small momentum transfers provided 
in [135H137] show that it is possible to determine both sizes of the nuclear core and halo. 
For the first time, the relevant experiments [161 H38[ 1139] were performed at GSI Darmstadt 
with the help of the ionization chamber IKAR |14U] developed at PNPI Gatchina. 

A. Experiment 

Differential cross sections da/dt were measured in inverse kinematics at GSI Darmstadt 
for proton scattering on nuclei of 4 ' 6 ' 8 He and 6 ' 8 ' 9 > n Li isotopes |1381 1139] . The measurements 
were performed at the equivalent proton energy E p w 0.7 GeV, the range of the momentum 
transfer squared t being 0.002 < \t\ < 0.05 (GeV/c) 2 . A primary 18 beam from the heavy- 
ion synchrotron SIS at GSI Darmstadt was focused on an 8 g/cm 2 beryllium production 
target at the entrance of the fragment separator FRS. The helium and lithium ions, produced 
by fragmentation of 18 nuclei, were separated by the FRS according to their magnetic 
rigidity. The intensity of the secondary He and Li beams was about 10 3 s™ 1 with a duty 
cycle in the range of 25 - 50 %. 

A scheme of the layout of the experiments is shown in Fig. 9. The main component 
of the setup is the hydrogen-filled ionization chamber IKAR, which served as a gas target 
and simultaneously as a recoil proton detector. IKAR was operated at 10 bar pressure, 
which insured the effective H 2 thickness of about 3 x 10 22 protons/cm 2 . IKAR consists 
of six independent identical modules. Each module is an axial ionization time-projection 
chamber, which contains anode plates, a cathode plate, and a grid (see insert in Fig. 9), 
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FIG. 9: Layout of the experimental setup for small-angle proton elastic scattering on exotic nu- 
clei in inverse kinematics. IKAR - hydrogen-filled ionization chamber which serves as gas target 
and detector of recoil protons. PC1-PC4 - multi-wire proportional chambers, which measure the 
projectile scattering angle. S1-S3 and VETO - scintillation counters for triggering and for beam 
particles identification. The ALADIN magnet with a position scintillator wall behind allows to 
identify the scattered beam particle and to discriminate the breakup channels. 



all electrodes being arranged perpendicular to the beam direction. The signals from the 
electrodes provided the energy of the recoil proton (or its energy loss in case it leaves the 
active volume), the scattering angle of the recoil proton, and the coordinate along the 
chamber axis of the interaction point in the grid-cathode space. 

The recoil protons in IKAR were registered in coincidence with the scattered He or Li 
particles. The momentum transfer could be determined either from the measured energy of 
the recoil proton or from the value of the projectile scattering angle, which was measured by 
a tracking detector consisting of 2 pairs of two-dimensional multi-wire proportional chambers 
arranged upstream and downstream with respect to IKAR (see Fig. 9). A set of scintillation 
counters was used for triggering and identification of the beam particles via time-of-flight and 
dE/dx measurements, while a circular-aperture scintillator VETO selected the projectiles 
which entered IKAR within an area of 2 cm in diameter around the central axis. 

In the case of the experiment with Li isotopes |139j a magnetic-rigidity analysis of the 
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scattered particles was also performed with the help of a large-gap magnet (ALADIN) and 
a scintillator wall behind it, which allowed one to exclude a contribution from the break-up 
channels. The systematic uncertainty in the normalization of the measured cross sections 
was estimated to be about 3%, while the uncertainty in the t-scale calibration was about 
1.5 %. 




0.02 
-t, (GeV/c) 



0.03 

2 



FIG. 10: Measured differential cross sections for p 6 ' 8 ' 9 - 11 ^ elastic scattering (dots) versus the four- 
momentum transfer squared. Solid lines are the cross sections calculated within the Glauber theory 
using the GO parameterization for the matter density distributions with the fitted parameters. 



The differential cross sections da/dt measured for the case of Li isotopes is shown in Fig. 
10. At first glance all cross sections have similar behavior. A steep rise of the cross sections 
with \t\ decreasing at \t\ < 0.004 (GeV/c) 2 is due to the Coulomb scattering. At \t\ > 0.005 
(GeV/c) 2 , the cross sections decrease with \t\ approximately as exponents. However, if one 
divides the cross sections by exponents (see in Fig. 11), then it is seen that the dependence 
of da/dt on \t\ for p 6,8 ' 9 Li and for p n Li scattering is different. The shape of da/dt for the 
case of proton scattering on 6 > 8 ' 9 Li nuclei (at 0.005 < \t\ < 0.05 (GeV/c) 2 ) is indeed close 
to that of exponents, whereas the shape of da/dt for p u Li scattering deviates significantly 
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from the exponential one. As will be discussed in Section VI B such a shape of the cross 
section is an indication for the core + halo structure of the investigated nucleus. 




0.01 0.02 0.03 0.04 0.05 
-t, (GeV/c) 2 

FIG. 11: Differential cross sections, the same as in Fig. 10, divided by exponents as explained in 
the text. 
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B. Analysis of the data 



The differential cross sections for proton elastic scattering on the studied nuclei were 
calculated by the Glauber formula (Eq. 2.8) using phenomenological density distributions, 
each having two free parameters, which were determined by fitting the calculated cross 
sections to the experimental data. As we discussed in Section [TTJ the pN scattering amplitude 
can be described by the standard high-energy parameterization 

/pN(q) = (iA;/47r)cr pN (l - ie pN )exp(-/3 pN q 2 /2) , (6.1) 

where a pp and the total cross sections of pp and pn interaction, e pp and e pn are 

the ratios of the real to imaginary parts, and (3 PP and /3 pn are the slope parameters. The 
slope parameters /3 PP and /3 pn were evaluated from the experimental data and partial wave 
analyses for free pp and pn scattering. 

Four parameterizations for phenomenological nuclear density distributions were applied, 
labeled as SF (symmetrized Fermi), GH (Gaussian-halo), GG (Gaussian- Gaussian) and GO 
(Gaussian-oscillator). In the SF parameterization |141j . the free parameters are the "half 
density radius" R and the diffuseness parameter a. The GH parameterization [16] is defined 
by the formfactor 

S(t) = (1 + az 2 )exp(z) , 

where R m is the root-mean-square radius of the matter distribution and < a < 0.4 and 
z = tR^/6. The GH distribution becomes a Gaussian one when a = 0, whereas for a close 
to 0.4 this distribution has a pronounced halo component. 

While the SF and GH parameterizations do not set apart the neutron and proton distri- 
butions, the GG and GO parameterizations assume that the nuclei consist of core nucleons 
and valence nucleons with different spatial distributions. The core nucleon distribution is 
assumed to be a Gaussian one in both the GG and GO parameterizations. The valence nu- 
cleon density is described by a Gaussian or a \p shell harmonic oscillator type distribution 
within the GG or GO parameterizations, respectively. 

The free parameters in the GG and GO parameterizations are the rms radii R c and R v 
of the core and valence nucleon distributions. The explicit expressions for the SF, GH, GG, 
and GO parameterizations are given in [20]. The studied nuclei were considered to have one, 
two, or four valence nucleons. The cores in 6 He, 8 He, 6 Li, 8 Li, 9 Li, and n Li were presumed 
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to have the nucleoli composition and the spatial structure similar to that of the 4 He, 4 He, 
4 He, 7 Li, 7 Li, and 9 Li nuclei, respectively. 




r, fm -t, (GeV/c) 2 

FIG. 12: The sensitivity of the differential cross sections for small-angle proton elastic scattering 
to the nuclear size and the radial shape of the nuclear matter distribution. 

The sensitivity of the different cross sections to the nuclear size and radial shape of the 
nuclear matter distribution is demonstrated in Fig. 12, where results of calculations for p u Li 
scattering at 0.7 GeV are displayed. The differential cross sections da/dt for small-angle 
elastic scattering is shown in Fig. 12 (b) as a function of t. They were calculated for a 
Gaussian matter distribution p{r) with different matter radii R m (see Fig. 12 (a)). A strong 
correlation between the slope of the cross section and the radius R m is obvious. The lower 
part of Fig. 12 demonstrates the sensitivity of the calculated cross section to the radial 
shape of the nuclear matter distribution. Two different nuclear matter density distributions 
are assumed, one being a Gaussian with R m = R c = R v = 3.5 fm, the other being the sum of 
two Gaussians - one for the core nucleons (with R c = 2.5 fm) and the other for the valence 
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TABLE VII: Root-mean-square radii R c , R v and R m correspondingly of the core nucleon, valence 
(halo) nucleon and total matter distributions, deduced from the data on small- angle proton elastic 
scattering. 



Nucleus R c (fm) R v (fm) 


R m (fm) 


4 He 




1.49(3) 


6 He 


1.88(12) 2.97(26) 


2.30(7) 


8 He 


1.55(15) 3.08(10) 


2.45(7) 


6 Li 


2.10(15) 3.00(32) 


2.44(7) 


8 Li 


2.48(7) 2.58(48) 


2.50(6) 


9 Li 


2.20(6) 3.12(28) 


2.44(6) 


n Li 


2.52(2) 5.98(32) 3.42(11) 



(halo) nucleons (R v = 6.4 fm). Both density distributions (shown in Fig. 12 (c)) have the 
same nuclear matter radius R m but significantly different radial shapes. Fig. 12 (d) depicts 
in the logarithmic scale the calculated cross sections related to these densities. In order to 
see more clearly the sensitivity of the cross sections to the density shape, the calculated 
cross sections are normalized by the exponential functions Coxexp(Bot). Here t = t — t*, 
t* = —0.01 (GeV/c) 2 , while the quantities Co and Bo are the values of the differential cross 
sections and their slope parameters at t — t*, corresponding to the two density distributions 
involved in the calculations. It is seen that the shapes of the cross sections is significantly 
different for the two cases considered, the shape calculated assuming the core+halo structure 
of n Li being similar to that of the experimental cross section shown in Fig. 11. 

The results of the data analysis using phenomenological density distributions are pre- 



sented in Table VII In the cases of 4,6 ' 8 He and 6 ' 8 > 9 Li nuclei, good data fit was achieved 
for all four density parameterizations applied. However, in the case of n Li, only GG and 
GO parameterizations, which allow different distributions of the core and valence (halo) 
nucleons, permitted good data description. 

The obtained results can be interpreted as an indication on the core+halo (skin) structure 
of 6 He, 8 He, 6 Li, 8 Li, 9 Li and n Li. Within the quoted errors the matter radii of the 6 Li, 8 Li 
and 9 Li nuclei are identical. This means that 8 Li and, especially, 9 Li are more dense nuclei 
than 6 Li. The latter is being considered to have an a + d spatial structure. The determined 
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core and halo radii represent clear evidence of a neutron halo in n Li. Indeed, the deduced 
halo radius = R v = 5.98(32) fm is larger than the core radius R c = 2.52 fm by a factor 
of more than 2. The total matter radius R m of n Li is significantly larger than those of 
the lighter Li isotopes. This result being in agreement with the data on nucleus-nucleus 
interaction cross sections discussed before. 

Among the studied nuclei of He and Li isotopes, 8 He and n Li have the most developed 
halo-like (skin-like) structure. In Fig. 13 (a) and (b), the core and total matter density 
distributions derived for 8 He and 11 Li within the GG and GO parameterizations are com- 
pared. In the case of 8 He it is still unclear if the valence neutron distribution is halo or skin 
like. The observed extended valence neutron distribution in n Li at the nuclear periphery is 
decidedly an outstanding halo. 




r, fm r, fm 

FIG. 13: Nuclear core and total matter distributions in 8 He and 11 Li, deduced from the cross 
sections for elastic p 8 He and p u Li scattering with the help of GG and GO density distribution pa- 
rameterizations. The shaded areas represent the envelopes of the matter and core density variations 
within the model parameterizations used, superimposed by the statistical errors. 

In the analysis of the cross sections with phenomenological densities discussed above all 
nucleon correlation (except the centre-of-mass correlations) in the nuclear many body density 
distributions were neglected. The deduced matter radii of 6 He and 8 He were i? m ( 6 He) = 2.30 
(7) fm and i? m ( 8 He) = 2.45 (7) fm [16]. However, the analysis of the same data carried out 
later by Al-Khalili and Tostevin in |142j using theoretical density distributions, the obtained 
radii of the same nuclei are larger by approximately 0.2 fm. Al-Khalili and Tostevin believed 
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that larger values of R m than those in [16] were obtained because few-body correlations 
in 6 He and 8 He were treated properly and used density distributions were with correct 
asymptotic. Authors stressed that the calculated cross sections and the deduced radii are 
somewhat sensitive to the few-body correlations in the many-body density and to the density 
asymptotic at large radii. This subject was considered in detail in [20J. It was concluded 
in contrast to [142] that there is very weak sensitivity of the calculated cross sections (at 
small scattering angles) to the nuclear correlations and to the density asymptotic. On the 
other hand, the calculated matter radii R m depend significantly on the density asymptotic. 
Theoretical density distributions [142j decrease with the radius increasing at large distances 
from the nuclear centre slower than the phenomenological density distributions in [IB] . We 
believe that provided reasons explain the larger values of R m obtained in [142] . 

The nuclear density distributions in nuclei with low binding energies should have long 
density tails. Proton elastic scattering is sensitive to the nuclear spatial structure including 
the most part of the halo, however it is not sensitive to small density tails which contain only 
of the order of 1 percent (or even less) of the total matter. A contribution of such tails to 
the value of rms matter radii can be estimated theoretically. An analysis of the pHe and pLi 
scattering cross sections with phenomenological density distributions including density tails 
taken from theoretical considerations was performed in [201 [139J. The rms matter radii for 
6 He, 8 He, and x lLi were found to be i? m ( 6 He) = 2.45 (10) fm, i? m ( 8 He) = 2.53 (8) fm, and 
_R m ( n Li) = 3.71 (20) fm. The relatively large error in the obtained value of R m in the case 
of n Li is mainly due to an uncertainty in the size and slope of the density distribution tail. 
A later analysis (not published) of the p u Li scattering cross sections assuming a smaller 
contribution of the density tail yielded i? m ( 11 Li) = 3.60 (20) fm. 

Combining the matter radii R m obtained from the data on proton elastic scattering with 
the proton radii R p , determined from the nuclear charge radii measured in laser spectroscopy 



experiments (see Table VI), the neutron radii R n and the thickness of the neutron skin (halo) 



5 np = R n — R p can be determined. Values for R p , R Q and <5 np are provided in Table VIII 



First experiments on the p 6,8 He and p 8 ' 9 ' n Li scattering in inverse kinematics [T6J I138[ I139j 
have shown that the intermediate-energy small-angle proton scattering is a useful means of 
investigation of the matter density distribution in light exotic nuclei. Future measurements 
of the cross sections for proton elastic scattering at larger momentum transfers [143] will 
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TABLE VIII: Total matter radii R m including contributions from density tails, neutron radii R n 
deduced from matter radii R m and the charge radii R c h, and the neutron skin (halo) thickness <5 np 
in the 6,8 He and 8 > 9 ' n Li nuclei. 

Nucleus R m , fm R n , fm <5 np , fm 
6 He 2.45(10) 2.68(14) 0.76(14) 
8 He 2.53(8) 2.73(10) 0.92(11) 
8 Li 2.50(6) 2.68(9) 0.52(10) 
9 Li 2.44(6) 2.59(9) 0.48(9) 
n Li 3.60(20) 3.96(25) 1.58(25) 

provide more detail information on the internal spatial structure of the studied nuclei. 

VII. CONCLUSION AND OUTLOOK 

There is an evident interest to the properties of radioactive nucleus. In nature there 
are 283 stable or very long lived nucleus [144 J and about (6 -v- 8) x 10 3 radioactive nuclides 
[84, 144J. Only about half of the existing nuclei have been studied so far. The nuclear size is 
one of the basic parameter of the nuclear density distribution studying which is important 
for understanding of the nuclear properties. 

A lot of information on the matter radii has been obtained from the nucleus-nucleus 
interaction cross sections, which can be measured for very low intensity beams of exotic 
nuclei. However, the deduced matter radii are somewhat model dependent and are subject to 
some uncertainties appearing due to approximations used in the calculations of the reaction 
cross sections. The optical-limit approximation significantly overestimates the calculated 
reaction cross sections, especially in the case of halo nuclei. The rigid-target approximation 
provides somewhat more accurate results. In principle, the reaction cross sections can be 
expressed using the exact Glauber theory formulas and then numerically calculated using 
Monte-Carlo technique 

Provided the intensity of the nuclear beams is sufficient, the matter density distributions 
in exotic nuclei can be studied in intermediate-energy proton elastic scattering experiments 
in inverse kinematics. The charge radii of exotic nuclei are measured very precisely with the 
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laser-spectroscopy technique. 

New experimental facilities for studying the properties of nuclei far from stability are 
planned to be built in the near future in Europe, Japan and the USA |!45H153j . Studying 
nuclei at these facilities will significantly increase our understanding of the spatial structure 
of the nuclei far from stability. 

Ambitious project NuSTAR at FAIR (Darmstadt, Germany) [147H149] will provide new 
fascinating possibilities for studying radioactive nuclei. New facility will produce intensive 
intermediate-energy and low-energy beams of nuclei far from stability which will allow to 
carry out versatile investigations of the nuclear properties. In particular, matter radii will 
be determined for long isotopic chains of many elements from the measured nucleus-nucleus 
interaction cross sections. Experiments on proton elastic scattering will be used to obtain in- 
formation on matter density distributions. Collaboration also plan to use laser-spectroscopy 
technique to measure nuclear charge radii and electron elastic scattering experiments to 
determine nuclear charge distributions. 

We expect that the improvement of the existing experimental technique and new facility 
will advance our understanding of the unstable nuclei. 
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[155] Many of unstable nuclei, for example n Li [53] have no excitation levels. 

[156] Interaction cross section is defined [15] as the total cross section for the processes of nucleon 
(proton and/or neutron) a^B removal from the incident nucleus. The difference between 
a^ B and a^ B was estimated [Ml 09] to be less than a few percents for a beam energy higher 
than several hundred MeV per nucleon and this difference will be considered in details in 
Section [TV] 

[157] We neglect the difference of electromagnetic and strong interaction nucleon radii. 

[158] Later the interaction cross sections of neutron-rich nuclei were described by Glauber Theory 

calculations with nuclear densities obtained from Skyrme Hartree-Fock calculations with 

accounting for more detailed picture, see, for example, |66j . 
[159] When using these data, it is necessary to keep in mind the criticism presented in |81j . 
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